Nonlinear sliding friction of adsorbed overlayers on disordered substrates 
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We study the response of an adsorbed monolayer on a disordered substrate under a driving force 
using Brownian molecular-dynamics simulation. We find that the sharp longitudinal and transverse 
depinning transitions with hysteresis still persist in the presence of weak disorder. However, the 
transitions are smeared out in the strong disorder limit. The theoretical results here provide a 
natural explanation for the recent data for the depinning transition of Kr films on gold substrate. 



PACS numbers; 68.35.Gy, 68.35.Rh, 81.40.Pq 

I. INTRODUCTION 

Adsorbed overlayers under a driving force have at- 
tracted much interest i-ecently^'^'^'^''^-^^i2iSi2ii24ii as im- 
portant examples of driven lattice systemsiSiii. The 
nonlinear response of the overlayer in the presence of 
pinning centers resulted in various steady states, cor- 
responding to a variety of different nonequilibrium dy- 
namic phases. The nature of these dynamic phases 
and the transitions from one phase into each other are 
the microscopic factors controlling the macroscopic slid- 
ing friction and boundary lubricatiorii^iiSiiSiii. Other 
closely related systems include vortex lattices in type- 
II superconductorsi2ii2i2242ii22i22i, charge-density wave^ 
and colloidal crystal layers2^. In general, both a peri- 
odic potential and disorder due to quenched defects are 
present in theses systems. In particular, in the sliding 
friction and boundary lubrication of surfaces the under- 
lying crystal lattice provides the periodic pinning poten- 
tial which competes with defect pinning centers due to 
impurities. 

For periodic substrates, it has been shown through 
molecular dynamics simulation studies that there are dif- 
ferent dynamical phases for the driven overlayer such as 
pinned solid, sliding solid and liquid phases. The transi- 
tions between these phases often involve strong hystere- 
sis effects^i, leading to important macroscopic behavior 
as the stick and slip motion, which can be observed ex- 
perimentally with the surface force apparatusi^ii^. Most 
of the theoretical results were based on two-dimensional 
models of overlayers with Lennard-Jones potential inter- 
action between the adatoms«^. The results on these dy- 
namical phases have been found to be universal and are 
reproduced even in models with pure elastic interacting 
potential^. More recently^, it has also been shown that, 
in the sliding state, the response of the overlayer to an 
additional force in the transverse direction can also lead 
to a depinning transition at low temperatures. Different 
sliding states also occur in presence of quasiperiodic sub- 
strates and have been investigated in detail recently for 



one-dimensional model o^i^° . However, the effects of dis- 
order on the main features observed for two-dimensional 
systems, which can arise, for example, from point defects 
on a initially periodic substrate, have not been investi- 
gated satisfactorily so far. Considerable analytical and 
numerical results are only available for driven lattice sys- 
tems described by overdamped dynamicsiSii^. Unfortu- 
nately, it is known that inertial effects significantly mod- 
ify the nature of the depinning transition-. Hence the 
results for these overdamped systems cannot be directly 
applied to the sliding friction problem involving adsorbed 
overlyers on a substrate. 

In this work, we study numerically the nonlinear slid- 
ing friction of adsorbed overlayers on disordered sub- 
strates with a specific model of point defects, where the 
strength of the disorder can be varied by changing the 
concentration of the defects. We investigate the longitu- 
dinal and transverse response for weak and strong disor- 
der by Brownian molecular dynamics simulations. We 
find that the main features, such as hysteresis, longi- 
tudinal and transverse depinning transitions previously 
observed for periodic surfaces^*?:, still survive in the pres- 
ence of weak disorder but are smeared out in the strong 
disorder limit. Recently, the depinning transition of Kr 
films on gold has been studied using the quartz crystal 
microbalance technique—. Our results for the longitudi- 
nal response bear a strong qualitative resemblance to the 
experimental data of the depinning transition. In par- 
ticular, it provides concrete theoretical support to the 
speculation that the aging effect on the depinning tran- 
sition is due to the increasing disorder of the substrate. 



II. MODEL AND SIMULATION 

The model of interacting adatoms on a substrate with 
point disorder subject to an external driving force for 
the present study is a simple extension of the one previ- 
ously used for periodic substrate a^i'^i^ . The dynamics is 
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described by the Langevin equation 

dU dV , ^ 
mri+mrir, ^ h f » + F (1) 

where ri = (xi,yi) is the adsorbate position, U = 
u{ri) is the substrate pinning potential, V ~ 
Tlii^j u(|ri — rj) is the interactions between adatoms, F is 
the uniform external force acting on each adatom and 
is a stochastic force, with zero average, and variance re- 
lated to the microscopic friction parameter -q, the mass of 
the particles m and the temperature T by the fluctuation 
dissipation relation 

<fr{t)fT{t')>='^VmkT5{t-t'), (2) 

where a represents Cartesian components. We choose a 
substrate potential which provides a periodic pinning for 
the overlayer except near a pinning center, given by the 
form 

u{y) — Uo [2 — cos(27ra::/a) — cos(27r?//a)] (3) 

Near a pinning center (point defect) located at Rp = 
{Xp,Yp), the substrate potential is locally modified to 
the form 

u(r) — Up [2 — cos{2TTx/a) — cos(27ry/a)] (4) 

for 

Xp - a/2 < X < Xp + a/2 
Yp~a/2<y<Yp + a/2, (5) 

The strength of the random pinning potential Up is 
chosen to be uniformly distributed in a range of 4[/o with 
an average < Up >= Uo- The position of the pinning cen- 
ters Rp — (Xp, Yp) is randomly chosen on a square lattice 
of lattice parameter a and linear size L. With this choice, 
a pinning center consists of a square region with a force 
that vanishes smoothly at the boundary. For Up — Uo, 
a constant, the entire substrate potential reduces to the 
periodic cosine potential in Eq. |S1 studied previouslyi*^. 
The strength of disorder is controlled by varying the con- 
centration Xd of randomly chosen pinning centers. In the 
dilute limit oi Xd << 1, corresponding to weak disor- 
der, this pinning potential has a predominant Fourier 
component at the wavevector k = {2rmT/a,2n'K/a). In 
the other limit of dense pinning centers with a;^ ^ 1, 
it represents an amorphous substrate with strong disor- 
der. The interaction between adatoms is represented by 
a Lennard- Jones pair potential 

vin,)^e[iro/n,)'^-2irJn,f] (6) 

where e is the well depth and tq is the particle separation 
at the minima in the pair potential. We have chosen the 
parameters in the interaction potential to be Vo/a ^ 1.5, 
e = Uo, and the adsorbate coverage at d — 1/2. This 
choice is made such that the ground state of the over- 
layer in the periodic substrate {xd = ) corresponds 



to the commensurate pinned c(2 x 2) structure studied 
previously. For the periodic substrate, the longitudinal 
and transverse response of the commensurate overlayer 
to an applied force in the present model have been exten- 
sively studied in connection with the sliding friction and 
boundary lubrication^'^'^ problem. Thus, this model pro- 
vides a convenient starting point for the present study of 
disorder effects. In the discussions below we use dimen- 
sionless units normalizing the force by 2ttUo, velocities 
by 2TrUo/r] and temperature by Uq/Ub- We study the 
nonlinear response of the overlayer at different temper- 
atures and pinning center concentrations via Brownian 
molecular dynamics simulationsSS.. Systems consisting 
of = 50 to 800 adatoms were considered with the 
time variable discretized in units of 5t = 0.002 — O.Olr, 
where r = [ma^ /Uo)^^'^ ■ The main results were ob- 
tained for e — Uo and rj — 1. For each set of parame- 
ters, calculations were performed for typically 1 — 4 x 10^ 
time steps to allow the system to reach steady state. 
Then an equal number of time steps are used to eval- 
uate time-averaged values of various physical quantities 
in this steady state. Finally, the results are also averaged 
over a number of (typically 6-11) different realizations of 
the disorder through the distribution of the Up value for 
the pinning centers. 



III. RESULTS AND DISCUSSION 

A. Depinning Transition and Longitudinal 
Response 

The longitudinal response to an applied force at a tem- 
perature T = 0.2, well below the melting transition tem- 
perature of the c(2 X 2) phase on a periodic substrate, 
Tm = 1.2, is shown in Figs. ^ and [3 for dilute Xd — 0.1 
and dense disorder Xd — I, respectively. In the for- 
mer case, the main features are similar to that found 
for the periodic substrate^. The drift velocity Vx is es- 
sentially zero below a critical value of Fa (neglecting the 
tiny contribution from thermal creep motion which is al- 
ways present for nonzero force). When the driving force 
exceeds the critical value Fa, there is a sharp depinning 
transition where the overlayer starts to slide and the drift 
velocity increases rapidly. Eventually a sliding regime for 
Fx >> Fc is reached where the longitudinal sliding fric- 
tion, defined as fj = F^/Vx, is close to the microscopic 
friction, fj ~ rj. There exists a hysteresis behavior for 
the depinning transition with unequal critical forces Fa 
and Fh for increasing and decreasing external force Fx, 
corresponding to the static and kinetic friction forces, re- 
spectively. On the other hand, in the dense disorder limit 

= 1 of Fig. 13 the hysteresis is completely washed out 
leading to equal static and kinetic friction forces Fa = Fi,. 

Recently Carlin et aW- have studied the depinning 
transition under an external driving force for Kr films 
on gold. They employed the quartz crystal microbalance 
technique. In this technique, the quartz crystal is driven 
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FIG. 1: Drift velocity Vx (a) and effective temperature T* 
(b) of the overlayer as a function of the external force Fx, 
for dilute disorder Xd = 0.1. The results are for T = 0.2 
and system size L = 20. The direction of variation of Fx is 
indicated by arrows. 



FIG. 2: Drift velocity Vx (a) and effective temperature T* 
(b) of the overlayer as a function of the external force Fx , for 
dense disorder Xd ~ 1. The results are for T — 0.2 and system 
size L = 20. The direction of variation of Fx is indicated by 
arrows. 



to resonance by an external rf voltage with the Kr film 
adsorbed on the gold electrodes. The driving force on 
the admolecules is proportional to the oscillation ampli- 
tude A of the electrodes and the resultant response of the 
Kr film is measured through the " slipping time" r which 
is inversely proportional to the sliding friction fj. For a 
given coverage of the Kr film, they found a sharp depin- 
ning transition when the oscillation amplitude exceeds a 
critical value. This transition shows the same kind of 
hysteresis effects as obtained in our weak disorder result 
shown in Fig^ A very interesting result they found is 
that the dcpinning transition changes as a function of 
aging of the system. As shown in Fig. the dcpinning 
transition is sharp and shows hysteresis behavior after 
1 day. However, the corresponding results after 10 days 
show a rounded transition without hysteresis. The au- 
thors speculated that this change is probably due to the 
deposition of foreign molecules onto the gold electrode 
substrate, acting then as defect pinning centers. In other 
words, the disorder of the substrate increases with aging. 
In Fig. we combine the numerical results for Xd = 0.1 
and Xd = I from the present model to show the effect of 
increasing density of defect pinning centers on the dcpin- 
ning transition. Our results bear a strong qualitative re- 
semblance to the experimental data shown in Fig. In 



particular, while the weak disorder (x^ = 0.1) substrate 
shows a sharp dcpinning transition with strong hystere- 
sis effect, the dcpinning transition for the strong disorder 
{xd = 1) substrate is smeared out and there is no trace of 
hysteresis in this case. In fact, our result for the strong 
disorder case bears a strong qualitative resemblance to 
the experimental data aged after 10 days. Thus, the 
theoretical results obtained in our present model study 
clearly support the speculationil that the aging effect on 
the dcpinning transition observed in the quartz microbal- 
ance experiment is due to the increasing disorder of the 
substrate. Another interesting feature of the experimen- 
tal results in Fig. 3a is the monotonic decrease of the 
slip time with aging. It should be noted that this effect 
is not seen in Fig. 3b because the nonlinear sliding fric- 
tion fj is normalized by the unknown microscopic friction 
parameter rj, contained in the effective model of Eq. (1). 
This parameter is essentially the frictional damping of 
individual adatoms due to coupling to the substrate. As 
the morphology of the substrate changes, the frictional 
damping may change significantly. So the monotonic de- 
crease of the slip time in Fig. 3a can be accommodated 
within our model as resulting from an increase in this 
microscopic frictional damping with aging. The calcula- 
tion of this microscopic parameter will certainly require 
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FIG. 3: (a) Slip time r as a function of the amplitude of sub- 
strate oscillations, obtained from experiments on Kr/Gd films 
using the quartz crystal microbalance technique, by Carlin et 
al. — (b) Inverse of the nonlinear sliding friction f) (longi- 
tudinal), normalized by the microscopic friction parameter rj, 
as a function of the applied force for weak {xd = 0.1) and 
strongly disordered {xd = 1) substrate, obtained from the 
numerical simulations of this work. 



a detailed theoretical model for this mechanism includ- 
ing realistic potentials which are outside the scope of the 
present work. 



B. Nature of the sliding state 

To understand the nature of the sliding state at differ- 
ent sliding velocities, it is useful to introduce, just as in 
the periodic case^, the concept of an effective tempera- 
ture T* of the sliding overlayer. This is defined via the 
velocity fluctuation as 



m{< > - < V >^)/2. 



(7) 



For a meaningful definition of this effective temperature, 
one would expect that the velocity distribution should 
be Gaussian and the width of the distribution related to 
the effective temperature by Eq. |7| This is confirmed in 
Fig. ^ where it is shown that the velocity distribution at 
Fx > Fc is in fact well described by a Gaussian both for 
Xd = 0.1 and Xd = 1- 




FIG. 4: Velocity distribution at T = 0.2, = 1.5, for 
Xd ~ 0-1 (a) and Xd = 1 (b). Continuous line is a Gaussian 
fit. 



The behavior of the effective temperature T* under the 
driving force carries important information for the struc- 
ture properties of the sliding state. As can be seen from 
Fig. n in the dilute disorder case with Xd = 0.1, this ef- 
fective temperature is much higher than T in the region 
Fh < Fj. < Fc and decreases for F^ » F^, although 
never approaching T as in the periodic case. This sug- 
gests that for weak disorder, the region for Fi, < F^ < Fc 
is strongly disordered and corresponds to a dynamically 
melted liquid phase. For F^ » Fc, the shding state is 
a solid state with weak disorder. For the substrate with 
dense disorder Xd = 1, the effective temperature remains 
much higher than T even for large driving force, suggest- 
ing that in this case the overlayer remains in a liquid like 
or floating state even at large driving force and high drift 
velocities. It is interesting to point out that an expression 
relating T* —T and 77 — 77 proposed by Persson^i has been 
found to work well for a periodic pinning potential. Ac- 
cording to this relation, in the limit of large driving force 
Fx > Fc, fj should approach 77 and T* should approach 
T. Our results in Figs, ^andl^lshow that this expression 
breaks down in the presence of disorder. While fj does 
approach 77 in the large driving force limit, T* remains 
higher than T in this limit in presence of disorder. 

To support the above interpretation between the ef- 
fective temperature T* of the sliding overlayer and its 
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structure, we have evaluated the structure factor defined 
by 

Sik)^^[<Y,e^^<'^--'^-^>] (8) 

for large driving force at Fx — 1.5, T — 0.2 in both the 
dilute and dense disorder limit in the present model. The 
angular brackets denote thermal averages and the square 
brackets disorder averages. The results are shown in Fig. 
|S1 For a fully ordered configuration the normalized struc- 
ture factor S{k)/N = 1 at the peak for a wavevector cor- 
responding to the ordered structure. As expected, for 
both Xd = 0.1 and = 1 cases, this value is much 
smaller at the nonzero wavevectors confirming that the 
sliding layer in the presence of the disorder does not have 
full long range order. For Xd = 0.1 the structure factor 
shows clear prominent peaks at wavevectors correspond- 
ing to a commensurate c(2 x 2) structure, {k — (l,l)7r/a), 
consistent with a weak disordered sliding solid state. For 
= 1, corresponding to the strong disorder limit, one 
can identify in the structure factor still well defined, al- 
though weaker, peaks at the positions corresponding to 
an incommensurate hexagonal structure. However, the 
information about the initial periodic substrate no longer 
exists. The structure of the overlayer in this strong dis- 
order limit can be characterized as an incommensurate 
floating state both in x and y directions rather than as a 
liquid state. 

Finally, we like to examine how the results presented 
here depends on the details of the modelling of the dis- 
order and defect pinning centers. In this regard, it is 
interesting to compare the results presented here with a 
previous work by Persson^ where the disorder was im- 
plemented by fixing one of the adsorbate atoms. It was 
found that in this case, corresponding to Xd = 0.0069 (lin- 
ear size L = 12), the effective temperature T* increases 
with increasing force for F > Fa and that the sliding in- 
commensurate solid phase observed for the periodic sub- 
strate is destroyed and replaced by a fluid state. At first 
sight this seems to be in total contrast to the results 
found in this work for the corresponding dilute disorder 
limit. To understand this difference, we note that in a 
strict two-dimensional model, all the adsorbate particles 
are confined to move in the x-y plane only. Hence, the 
structure of the sliding layer has to undergo strong dis- 
tortions in the vicinity of the fixed impurity atom, due 
to the hard core repulsion of the Lennard- Jones inter- 
acting potential, leading to a fiuid phase. Thus, in this 
fixed-impurity model, even a low concentration of impuri- 
ties actually correspond to strong disorder, in sharp con- 
trast to the model studied here. One would expect that 
even for the fixed impurity model, a very small concen- 
tration of impurities should lead to a more ordered sliding 
state. To verify this, we have computed the structure fac- 
tor for the fixed impurity model at two concentration of 
Xd — 0.0069 ( linear size L = 12) and Xd = 0.0012 ( linear 
size L — 30) in the steady sliding state at F.^ — 1.5. For 
L = 12 as shown in Fig. |^, the corresponding structure 




FIG. 5: Structure factor S{k)/N at F^, = 1.5, Fy = for 
dilute disorder a;^ = 0.1 (a) and dense disorder Xd ~ 1 (b). 
The results are for T = 0.2 and system size L = 30. 

factor is characteristic of a fluid phase, in agreement with 
the findings of Persson"^. However, for the larger system 
size L = 30 in Fig. (Sja, and therefore more dilute disor- 
der, sharp peaks develop corresponding to an hexagonal 
structure. This structure is very similar to what we found 
in the present model for = 1 as shown in Fig. \Sjp in- 
dicating that the overlayer in this situation corresponds 
to an incommensurate hexagonal floating phase. 



C. Transverse response at large driving force 

In this section, we focus our attention on the nature 
of the sliding state at large driving forces F^ > Fc, and 
its relation to the response of the sliding layer to an ad- 
ditional transverse driving force. It turns out that the 
nonlinear response of the overlayer in the transverse di- 
rection is sensitive to the nature of the sliding state, thus 
serving as an ideal probe for the structure of the sliding 
state which is otherwise difficult to study. 
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FIG. 7: Transverse velocity Vy (a) and effective temperature 
T* (b) of the overlayer as a function of an additional force 
along the y-axis at fixed F^, for dilute disorder Xd ~ 0.1. 
Arrows indicate the direction of force variation. Results are 
for Fj; = 1.5, T = 0.2, and system size L = 20. 



FIG. 6: Structure factor S{k) /iV at T = 0.2, = 1.5, = 
for the periodic substrate and fixing one of the adsorbates at 
the potential minima for a system size L = 12 (a) and L = 30 
(b). 



We have studied the transverse response of the over- 
layer in a steady sliding state along the x direction. For 
a fixed large driving force = 1.5 > Fc, an additional 
variable transverse force Fy is applied, and the average 
transverse drift velocity Vy as a function of Fy is deter- 
mined. The resultant transverse response is shown in 
Fig. [7|and|Sl. For dilute disorder Xd = 0.1, we find that 
there is a transverse depinning transition with a critical 
transverse force Fya- This transverse depinning transi- 
tion is accompanied by hysteresis behavior with the lower 
critical force Fyt, ^ 0. The existence of this transverse 
depinning transition for weak disorder is consistent with 
our finding of the nature of the sliding state at large 
longitudinal driving force. As shown in Fig. ^ for the 
shding state with Fy = and F^ > F^ we have T* ~ T. 
The structure factor shown for this state in Fig. Ohas 
prominent peaks at the c{2x2) commensurate structure 
despite the disorder resulting from a random distribution 
of the pinning centers. Thus, even while sliding along the 



X direction, the overlayer can be pinned in the transverse 
direction by the commensurate pinning potential. This is 
origin for the existence of the transverse depinning tran- 
sition. In the other limit of dense pinning centers Xd — 1 
corresponding to strong disorder, the transverse response 
in Fig. |S1 shows no hysteresis and the transverse crit- 
ical force Fya appears to vanish. Again, this behavior 
can be understood from the nature of the corresponding 
sliding state at Fy = 0. As shown in Fig. |21 in this 
regime of strong disorder, T* » T when F^ > Fc- Here 
the overlayer film is fully disordered and incommensurate 
with the substrate along both the x and y directions. As 
a result, there is no pinning force along the transverse 
direction, hence the absence of a transverse depinning 
transition. 



IV. CONCLUSIONS 

We have studied the nonlinear sliding friction of ad- 
sorbed overlayers on disordered substrates with point de- 
fects (pinning centers). A simple model is introduced 
which is a generalization of the driven monolayer model 
previously used for periodic surfaces-*'^'-''. The present 
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point defects on an initially ordered substrate. We have 
studied the longitudinal and transverse response in this 
model by Brownian molecular dynamics simulations in 
both the weak and strong disorder limit. The results 
show that the main features, such as hysteresis, longi- 
tudinal and transverse depinning previously observed for 
periodic surfaces^il, still survive in presence of weak dis- 
order but are washed out in the strong disorder limit. 
The results for longitudinal response bear a strong quali- 
tative resemblance to the recent experimental data of the 
depinning transition of Kr films on Gd, using the quartz 
crystal microbalance techniqueii. This qualitative agree- 
ment supports the interpretationii of the experimental 
data that the aging effects on the depinning transition 
are due to the increasing disorder of the substrate. The 
present model can also be used to gain insight into other 
interesting features of the experimental results such as 
the coverage dependence of the depinning transition and 
sliding friction. We expect that the coverage dependence 
in the weak disorder limit to be similar to that obtained 
in numerical simulations for a periodic substratei but 
quite different in the strong disorder limit. Moreover, 
within the present model, a study of the nucleation of 
topological defects and detailed microscopic mechanism 
of the depinning transition itself should be infesting for 
driven lattice systems in general. 



FIG. 8: Transverse velocity Vy (a) and effective temperature 
T* (b) of the overlayer as a function of an additional force 
along the j/-axis at fixed F^, for dense disorder Xd ~ 1- Arrows 
indicate the direction of force variation. Results are for Fx — 
1.5, T = 0.2, and system size L = 20. 



model incorporates in a convenient way the effects of 
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